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Abstract 

A quasiclassical correspondent for the fermion degrees of freedom is obtained by using a time- 
dependent variational principle with Grassmann coherent states as trial functions. In the real 
parametrization provided by the canonical coordinates, these states satisfy a closure relationship, 
and this is used to calculate the partition function at finite temperature. The particular exam- 
ple considered here consists of a many-fermion system in a quantum double-well. Entanglement 
between the one-particle orbital states leads to deviations from the normal Fermi-Dirac distribu- 
tion. This deviation is reflected in the dependence of the chemical potential on concentration. In 
the physical example of two adjacent CuO planes in a high- Tc superconductor, the dependence is 
suggestive of the pseudogap observed at temperatures T > Tc . 



1 



I. INTRODUCTION 



The classical limit of an elementary quantum dynamical system can show spectacular 
variation in the nature and number of the degrees of freedom. While the quantum dynamics 
is mainly stationary, or in transitory regimes given by a linear superposition of an infinite 
number of harmonic-oscillator modes (the stationary states), the classical dynamics is in 
general nonstationary and nonlinear, and it involves a relatively small number of degrees of 
freedom. 

The transition from quantum to classical statistics in a many-particle system at finite 
temperature, resulting from a change in the underlying microscopic dynamics from quantum 
to classical, is less well explored. The usual approach is based on the Wigner-Kirkwood 
expansion of the quantum partition function in powers of the Planck constant h. However, 
convergence of this series is assured only when the thermal energy fc^T is large with respect 
to the level spacing 

Of course, the classical Maxwell-Boltzmann statistics can be recovered from the quan- 
tum Bose-Einstein or Fermi-Dirac statistics in the limit of high temperatures. However, 
this limit is not exactly a classical limit because the constituent particles remain quantum 
objects. Moreover for spin-i particles, a classical correspondent exists only in the sense of 
a pseudomechanics, defined in terms of Grassmann variables 0. 

Classical degrees of freedom at the quantum level can be introduced with the concepts of 
a molecular (or nuclear) mean field or of spontaneous symmetry breaking. In such situations, 
the infinite degeneracy of a metastable ground state is lifted by the emergence of a small 
number of collective (mean-field) degrees of freedom. At low energy, the collective motion 
is quantized, but with increasing energy, a transition to classical behavior may occur. 

The thermodynamics of quantum many-body systems can be affected not only by the 
classical behavior of certain degrees of freedom of the constituent particles, but also by 
macroscopic quantum coherence phenomena. Such phenomena occur, for example, when a 
quantum system is split by a potential barrier into two entangled components localized in 
distinct classical regions. This situation appears in the generic case of two quasi-2D (two- 
dimensional) layers of electron gas separated by a thin insulator. Examples of such bilayer 
structures that are particularly important for applications are provided by the CuO planes 
of high-Tc superconductors and by the tunnel junction of the recently developed 2D-2D 
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transistor ^. 

A suitable framework for understanding the quasiclassical behavior of quantum systems is 
provided by the time-dependent variational principle 0. This principle leads to quasiclassi- 
cal dynamics, obtained by constraining the quantum dynamics from the infinite-dimensional 
Hilbert space to a finite-dimensional manifold of trial functions. The usual mean-field ap- 
proximations of many-body theory can be recovered in this framework by choosing trial 
functions that are coherent states of a semisimple Lie group. However, by the choice of the 
generators, these approximations provide the classical analog for particle-hole or particle- 
particle operators rather than for individual fermion fields. 

A direct quasiclassical correspondent of the Fermi operators can be obtained by using 
as trial functions the Grassmann coherent states 0, . In Sect. II it will be shown 
that by a suitable definition of Hermitian conjugation and of the product between conjugate 
variables, the Grassmann coherent states admit a real parameterization. With respect to this 
parametrization, the variational principle associates a classical Hamiltonian system for each 
fermion degree of freedom. The integral of the projector on the Grassmann coherent states, 
taken over the invariant volume element associated with the canonical variables, provides a 
resolution of the identity in the many- fermion Hilbert space (see Appendix). This important 
relationship was obtained previously ^ by direct integration over Grassmann variables. 
However, the approach based on the Berezin integral may lead to ambiguities, especially 
in the integration of multiple products involving dual Grassmann variables. As we show, 
the present definition removes the difficulties encountered in the semiclassical quantization 
of the phase integrals 0, [l^ and leads to the correct Bohr-Wilson-Sommerfeld formula. 

The resolution of the identity allows us to reduce the calculus of the trace in many-fermion 
Hilbert space to an integral over Grassmann coherent states, and this is particularly suited 
for statistical calculations. At the end of Sect. II, we show that the present approach leads 
to the standard Fermi-Dirac distribution function. 

The effect of quantum entanglement on the particle distribution functions in coupled 
Fermi systems is discussed in Sect. HI. The particular coupling selected here is tunneling 
through a potential barrier, but this is only a heuristic model. Similar results can be ob- 
tained by using any interaction that splits the degenerate levels. It is seen that the level 
splitting and the difference between the chemical potentials of the two components combine 
in a single quantity that determines the entangled distribution function. In the physical 
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example of the CuO planes in the high-Tc compound La2-a;Sra;Cu04, we show in Sect. IV 
that the entanglement brings the normal Fermi-Dirac distribution closer to one of a super- 
conductor. An observable directly related to the distribution function is the shift of the 
chemical potential as a function of concentration. It is shown that despite the simplicity 
of the model, the shift calculated for an entangled Fermi-Dirac distribution is in reason- 
able agreement with available experimental data. A summary of the main results, and the 
concluding remarks are presented in Sect. V. 

II. GRASSMANN COHERENT STATES AND PARTITION FUNCTIONS 

Consider a model space containing a set of L orthonormal one-particle fermion states 
{V'j, i = 1, . . . , L}, and let a\ (oj = j ) be the associated fermion creation (annihilation) 
operators in Fock space with respect to the vacuum state |0), Oj |0) = 0, such that ipi{x) = 
{x\ a\ |0), {ai,aj} = aiaj + a^aj = and |a|,aj| = 5ij. Coherent fermion states |C) can be 
defined as eigenstates of the fermion annihilation operator Oj, with eigenvalues Q [0, 

a^ 10 = 10 • (1) 
This definition, together with the anticommutivity of the fermion operators, leads to 

{a„a,}|O = ±{0,0}IO=0, (2) 

with the upper (+) sign on the right-hand side if and commute, [oj, 0] = 0, and 
the lower (— ) sign if they anticommute, {a-iXj} = 0. In both cases {QXj} = 0) which 
means that the "eigenvalues" Q, i = 1, . . . , L cannot be scalars but should be considered 
Grassmann variables. Thus, coherent fermion states exist only in a Grassmannian extension 
of Fock space. 

Let be the number of independent variables among Q, i = 1, . . . , L, and let {wa}, 
a = 1, . . . , N < L, be a basis in the complex linear space Fjy of the Q. The product of 
k < N distinct basis elements is an element of an (^)-dimensional linear space Q^. For 
example, Q% denotes the space of complex scalars and Qj^ = Tn . The direct sum 0^o^Af 
is the 2^-dimensional linear space of the full Grassmann algebra Qn- 

To calculate matrix elements of many-body fermion operators between Grassmann coher- 
ent states, we must define the Hermitian conjugation operation "|" not only for operators 
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but also for Grassmann variables. The operation is introduced here by using the linear 
space dual to Fat, that is the linear space of complex linear functionals on Fat. The basis 
in F^ dual to the basis {wa} of F^r comprises functionals defined in terms of unique 
scalar- valued mappings "■" such that ■ wp = wp ■ w^^ = 6a/3 JT^- This relation allows us 



to introduce the involution operator * : F^v — > F^, related to the Hermitian conjugation by 
iQa,y = ale;, such that for any e Fjv, (C*)* = 0, (00)* = OC, and (zQ)* = z*q, 
where z is any complex number and z* is its complex conjugate. 

It is often convenient to think of the scalar-valued mappings denoted by "■" as products 
in the direct-sum space Cn = F^ © Fat. Duality then assigns a complex scalar value to 
any product C* ■ = J2a=i Grassmann variables C,* = J2a=i ^f*"^! ^ and 



(k = Sq=i z^Wa G Fjy. The space generates a 2^^-dimensional Clifford algebra Ci 



p!3| , |T^, |T5[ in which the scalar product of null vectors and can be identified with the 
anticommutator: C* ' C = | {C*; 0} of their associative Clifford product. Physical quantities 
are evaluated as the scalar part of expressions in C£n,n, and unless otherwise noted, this 
restriction is assumed in expressions of expectation values below. 

When L = 1, an unnormalized Grassmann coherent state can be expressed as 

C) = e^<^< |0) = (l ± Cial) |0) = |0) ± Ci |1) , (3) 

where |1) = a\ |0) . Since the (scalar part of the) overlap 

|CIC) = 1 + 0-0 (4) 

is a real number, the normalized state (|l|) is 

(i±cA) |o) 

10 = \ ^ . (5) 

Vi + O-Ci 

It may be verified that since aia\ |0) = < [ |0) = |0) and (f = 0, 

Ciaial |0) Ci |0) 



ailO 



yr+CTO VI + ■ Ci 

Ci (i±Ci«l) |0) 



Ci 10 • (6) 



VI + ■ Ci 

In the L-level case it is convenient to assume that Q and a] anti-commute, and we therefore 
choose the "— " sign in Eq. (^. Then the normalized Grassmann coherent state has the 



general expression 

The expectation value of the particle number operator a|.afc in this state is 

= (CI aia, |C) = (CI C.^C. 10 = T^T^ • (») 

Note that C^Cfc is not a pure scalar in CC.n,n (it also contains a bivector part), and since 
the coherent state |C) is also not a Grassmann scalar, we do not restrict CfcCfc? but only 
its expectation value (C| C^Cfc |C) > to its scalar part. When = L, it is possible to choose 
Cfc = ZkWk- In this case, Q ' (k = ^l^k and the absolute value of the complex variable 
-2fc = Pk^^'^'' is related to the particle number by = p^/ (1 + p|)- This relation limits 
to the range from (when = ) to 1 (when p^ ^ oo ). With such parametrization, 
Grassmann coherent states become su (2) coherent states (see Appendix). The values of Uk 
are quantized by the anticommutation relation |a|, aj | = Sij or by the old quantum theory 
(see below). 

A quasiclassical approximation of the quantum many-fermion dynamics can be obtained 
by using the Grassmann coherent state |C) as trial function in the time-dependent variational 
principle. If |C) depends on time only through the complex scalar variables z^, then the best 
approximation for an exact solution of the time- dependent Schrodinger equation satisfies 

4 [\c\ihdt-H\Odt = o , (9) 

where H is the Hamiltonian of the many-fermion system, ti,t2 are fixed times, and 6x 
denotes the variation of the integral with respect to the real functions of time Pk{t) and 
(pk(t), k = 1, . . . ,L, denoted together by Xk, k = 1, . . . , 2L. 

The equations of motion determined by Eq. (^ have the general form 

1=1 

where H = (Cl-f^lC); and Uij = —ujji = 2h'^{{di( \dj()) are the coefficients of a nonde- 
generate 2-form (a symplectic structure) on the manifold of Grassmann coherent states. A 
canonical transformation to the new variables 

qk= \ j^^sm(f)k , Pk = \ cos (l)k , (11) 

V 1 + Pfc V 1 + Pfc 
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defined on the disk of radius ylh in the 'pq plane, reduces further the equations of motion 
to 

gfe = ^ , Pfe = ■ (12) 

If the states ■i/'j, i = 1, . . . , L are one-particle energy eigenstates with eigenvalues ej, then a 
many-fermion Hamiltonian of the form Hf — Yl!k=i ^kO^k^Lk can be considered as the Hamil- 
tonian of a "Fermi oscillator." Since 

Pk 

for the Fermi oscillator Hf = J2k=i ^k [1 — (pI + qI) and the equations of motion have 
the solution Pk{t) — const, (pkit) — —tek/h. This represents a harmonic oscillation of 
the variables pk and qk at the angular frequency Ck/h. Moreover, in this particular case, 
the Grassmann coherent state \0 (t) is not an approximation but an exact solution of the 
time-dependent Schrodinger equation, 

ih^ = Hf 10 . (14) 

The "classical" variables pk and qk, have been introduced without reference to the physical 
degrees of freedom of the particles (space coordinates -|- spin). However, it is interesting to 
note that they behave in all respects as real classical variables. In the pq plane the orbits are 
circles of radius \/2h (1 — n^). The Bohr-Wilson-Sommerfeld quantization formula selects 
only the orbits having an area which is an integer multiple of the Planck constant h, such 
that 

27rh (1 — Uk) — nh , n — 0,1, .... (15) 

This shows that 1 — rik should be a non-negative integer, and the only possible values of Uk 
are and 1. 

The remarkable property of Grassmann coherent states of providing a resolution of the 
identity operator in many-body Hilbert space (see Appendix) 

d'Q\C){C\ (16) 

makes them particularly suitable for applications in statistical mechanics. The parti- 
tion function Z of a quantum A^p-body system is defined by the trace in Fock space. 



Z = Tr{e~^^) = (H I e"''^ I [n]) , where P = l/ksT, H is the many-body Hamil- 
tonian operator, and \[n\) = \nin2...) represents one of a complete set of states, labeled 
by the occupation numbers of the single-particle levels, and restricted by the condition 
Ylk = . If these states are eigenstates of H with eigenvalues then Z = ^^^j e~^^["i . 
Moreover, if Efn] = X^fc^fe^fc, then Z = J2[n]Y[k^~^^''''^'' ■ ^he sum S[„,] can be calculated 
if the condition of having a fixed number of particles is relaxed and the are indepen- 
dent ll6|. The average number of particles can still be fixed by introducing the chem- 
ical potential /i; then Z = Ylk becomes the product of one-level partition functions 
Zk = Srifc exp [—Puk (cfc — /^)] with = 0, 1 for Fermi-Dirac statistics and Uk = 0, oo for 
Bose-Einstein. 

The statistical significance of the "Fermi oscillator" appears when the corresponding 
partition function Zf is calculated as an integral over the manifold of Grassmann coherent 
states. The resolution of the identity presented in the Appendix allows us to write 

L 



Z^ = Tr 



-p{Hf^,Nr)j = JJZfc (17) 
k=l 

where Nt = Ylk ^\^k is the total particle number operator and 

Zk = lj dn^ {Ck\ e-^'^I-'^^-^^'^-'^) \Ck) . (18) 

Because rffi^ = d^kdcj)^ = dp^dq^/h, and for any integer power n, ^a^Ofej \(k) = 
(Cfcl olak \Ck) = rik with nk = l- {pl + ql) /2h, we get 

= 1 + e-'3(^*-'^) (19) 

recovering the standard result for fermions. 

If efc — /i, k = 1, . . . ,L, and T are such that \ek — fi\ « ksT, then e~^(^'=~^^ — 1 ^ 
—(3 (e/c — fi). In this case, Zf takes the form of a classical partition function, defined by the 
phase-space integral of exp [— /3 (Hf — fJ^Af)] with respect to the volume element d^Q, where 
Af is the expectation value of Nj- . 
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III. DISTRIBUTION FUNCTIONS IN ENTANGLED FERMI SYSTEMS 



The interplay among classical localization, quantum coherence, and temperature on the 
one-particle distribution function plays a central role in the many-fermion fusion or fission 
reactions and in the operation of tunnel junctions. Localization occurs when a container 
filled with Fermi gas is split into two identical, classically distinct components. The opposite 
process, when two separated, identical Fermi systems merge to form a single large system, 
represents delocalization. In both situations, the intermediate stages can be modeled by two 
entangled Fermi systems confined in classically distinct potential wells that are separated 
by a finite barrier. 

The problem of a system of fermions placed in a double-well potential at zero temperature 
was studied previously in the context of the predicted exchange-driven spontaneous transi- 
tions from bilayer to monolayer ground-state configurations Jl^. The influence of coupling 



on the distribution function at finite temperature is studied here by considering a system of 
fermions of mass M placed in a one-dimensional potential, consisting of two identical rect- 
angular potential wells of width w and depth [/, separated by a distance d. The bounded 
one-particle energy levels in this case are = xlU, where x = Xk satisfies 



tan (wx) = 2aVl x 

with w = w ^2MU/h? and d = d^2MU/h?. For the fully merged system {d = 0, = 2w) 
as well as for each single well {d = oo, Wd = w), the solutions of 



tan [wdx) = - 1 ' ^ ^ 

Let us denote by \kl) and \k2) the orbital wave functions localized in wells 1 and 2, respec- 
tively, which at (i = oo correspond to the same energy level e° = xlU. At large but finite d, 
tunneling hfts the degeneracy. The eigenstates of the coupled system are represented by the 
entangled symmetric and antisymmetric superpositions \k±) = {\kl) ± \k2))/\/2, separated 
in energy by the tunnel splitting = — e^^ ^ 2V^(/c), where 

V,{k) = ^e-'^'v^ (22) 

is the tunneling matrix element. The Hamiltonian of a many-fermion system distributed 
over the levels of the entangled wells is ife = '^km ^km, where m = ±1/2 is the z-component 



of the spin, and 

hkm = Cfe aLm"fe-"» + 4'^l+m(^k+m (23) 

/ 1 t \ ( ] t ^ 

= [^Ikm'^^km + a2km(^2km j + ^ (^Olfem'^Sfcm + ^fem'^l^m j " 

In a quantum description of the entangled system, the corresponding partition function is 
Ze = T^(^e~^^^''~^^'^^). However, as long as the two wells remain classically distinct near 
equilibrium, we can assume that each one develops its own chemical potential such that 
Zf. = Tr(e^'^*^-^'=^'^^^i~^2^2)). Following the procedure presented in Sect. II, Z^, can be 
reduced to an integral over Grassmann coherent states, 

Z, = ^ j d^l^id^Oa (CI e-^(^«-'^i^i-'^2^^) IC) 



"2 '"2 

k=l 



where 



and 



exp (-CunAlrn " C2A„/4a„,) 

ICfcm) = \ ^ ^ |0) . (26) 

V (1 + Clfem ■ Clfcm)(l + C2fem ' C2fem) 



The operator exp 



has the expansion 



1 + (e"'^ COsh% - l){a\^^aihm + 4fem«2fem) 

rjk 

+ Z^(«lfcm«2fcm + 4fem«lfem)] (27) 
'Ik 

+ (e^"* - 2e'''= cosh77fe + l)a\f,^aikm(Akm'^2km , 



where 0;^ = -/3(e^ - (//i + /i2)/2), rjk = a/sT+t^, = /3Afe/2, and r = /3(//i - //2)/2. 
Denoting by ((A)) the integral 

((A)) = y dflikmdn2km iCkm \A\ (km) (28) 
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we get 



((1)) = 4 , = {{al,^a2km)) = 2 , (29) 

such that Zkm has the expression 

= 6^"'= + 2e"'^ coshr^fc + 1 . (30) 

The quantity that contains the effect of entanglement on the distribution function is ?7a; = 
\/sl + r^. If \sk\ « \r\, then Zf^^n becomes the product ZkimZk2m of the Fermi distributions 
for the orbital levels \kl) and \k2), 

Zkim = l + e-^^^°-^^^ (31) 
Zk2m = l + e-''^'°-'''\ 

When |r| << \sk\, then Z/^m is the product Zk_mZk+m of the distributions 

Zk_m = 1 + e-^^^^-f^^ (32) 

and fi = (yUi + /i2)/2 is the Fermi level of the fully merged system. 

For the electron layers in the 2D-2D tunneling transistor, the difference between the 
chemical potentials is related to the applied voltage V, /ii — /i2 = eV . In the case of two 
adjacent CuO layers of a high- Tc compound, a small difference may appear due to the slight 
asymmetry observed in the placement of the oxygen acceptors between the planes [|T8|. 



IV. ENTANGLED FERMI DISTRIBUTIONS IN HIGH-Tc SUPERCONDUC- 
TORS 

In high-Tc materials the charge carriers are confined to the CuO planes, but the interlayer 
coupling is important in both the normal and superconducting states. At intermediate tem- 
peratures and low doping, the interlayer coupling is dominated by the magnetic interaction 
of the electron spins, producing antiferromagnetic ordering. The charge transport between 
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the layers, practically absent in this state, is assumed to have both a coherent, momentum- 
conserving component and an incoherent part due to impurity-assisted hopping |jl9[. With 
increasing doping, at lower temperature the coupling takes the form of coherent Josephson 
pair tunneling, supposedly one of the basic phenomena responsible for the superconductivity 
at high temperatures. In this perspective, the observed decrease of transition temperature 
Tc with an increase of the number of CuO planes seems to be explained by the loss of the 
phase coherence [|20|. 



Layered copper oxide compounds are peculiar not only in their high values compared 
to conventional superconductors, but also in various anomalous properties observed in their 
normal state. Anomalies showing a common trend with respect to the temperature and to 
concentration of charge carriers are confirmed by a wide variety of experimental techniques 
||21|| . Many anomalous features can be interpreted in terms of a "normal-state gap" or 
pseudogap [^, ^ ^ of unknown origin that turns into the BCS gap when the material 
becomes superconducting. 

The results presented in Sect. Ill indicate that entanglement between adjacent layers 
can produce significant deviations from the standard Fermi-Dirac distribution that may 
contribute to the observed anomalies. The magnitude of such effects is evaluated here within 
a schematic model in which the energy of the electrons in a CuO plane without coupling 
has the simple form 

4 = e° = 6, + 6^ , (33) 

where k = (e^, n) denotes the couple of a continuous index and a discrete index n. These 
indices correspond to the terms ep and e'^ of the total energy due to the motion of the electrons 
along directions parallel and normal to the planes, respectively. The in-plane energy term 
is approximated by a single band, G [emin,emax], with -e™„ = emax = 1 eV, and a 
Fermi level ep = —0.45 eV [^. Confinement is modeled by a one-dimensional potential 



well placed normal to the planes. As suggested by the expression for the hole binding 
potential in oxygenated La2Cu04+y, |T^, the potential well is about 35 meV deep and has 
a radius a = 8 A. Thus, is calculated by using the square- well estimate = —(1 — x^)U , 
Xn ^ nrc /w < 1, where f/ = 35 meV, and w = 2a = IQ A. The interlayer coupling matrix 
element tc depends on doping, and in La2-xSr^.Cu04, tc ~ O.GGx'^'^'^^ eV |]2^. This shows 



that tc depends on the hole concentration 5 ~ a; and on the band filling / = 1 — 5. 
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The anomalous behavior attributed to the formation of a pseudogap occurs in general in 
underdoped compounds. In La2-xSr2;Cu04, Tc attains a maximum of ~ 40 K at Xopt ~ 0.15. 



Measurements of the specific heat coefficient |]25| in this compound indicate the existence 



of a pseudogap for x < 0.22. In the range 0.1 < s < 0.22, the pseudogap has an almost 
linear decrease with x, from ~ 200 meV to ||2B|. At x ~ 0.1, when the pseudogap is 



large, the interlayer coupling strength is tc ~ 1 meV. This value is large compared to the 
one obtained from Eq. (^2|) , indicating that it may include contributions due to other inter- 
actions, stronger than tunneling coupling. In the present calculations, the energy splitting 
Ak = An = 2Vc(n) was calculated by taking Vc{n) = 0.001 + 1.5^/x^ eV. This includes the 
realistic term tc = 1 meV, and a term weakly increasing with the excitation energy along 
the normal to the plane. 

The number of particles A^^i in the CuO plane 1 is related to the partition function Z^m 
of Eq. dag) by 

Here is calculated as an integral over the planar band and a sum over the bounded 
energy levels for motion normal to the plane, such that 

Y,^Y.^,S f de, , (35) 



where Np = m*/(27r^^) is the planar density of states, m* is the effective band mass p8 



and S is the plane surface. Thus, the planar density Dp = Ni/S of charge carriers takes the 
form 

Dp= / rfe(n)(,) , (36) 



where 

gOen _|_ 2 cosh rjn + e 



\ >r^e°^" + cosh r/„ + r/r/„ sinh r/„ 
n)/ N = 2Np > , (37) 



"en = -/5(e + e^) + /5(/ii + Ai2)/2, = A/sfTr^, s„ = /5A„/2, and r = /5(/ii - /i2)/2. 
The density is shown as a function of e at T = 120 K and r = in Fig. 1 (c). For 

comparison, on the same plot we display (a) the standard Fermi-Dirac distribution obtained 
in the absence of interlayer coupling (Vc = 0, s„ = 0), and (b) the distribution 

^ E,(e/^^« + i: 



n 
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of a superconducting system at finite temperature [^. In this expression, tlie quasiparticle 
energy = a/ (e° — ep^ + A^^ was calculated by using a fictitious superconducting gap 
Asc = 0.6 eV. 

The comparison between these functions shows that although different, as far as devia- 
tions from the normal Fermi-Dirac distribution are concerned, (b) and (c) may account for 
similar effects. An observable related directly to the particle distribution is represented by 
the shift AyU of the chemical potential as a function of concentration. This shift was sub- 
ject to a detailed experimental investigation, observing that in La2_xSr2.Cu04, an anomaly 
occurs at low doping levels [2^ . 

The comparison of the experimental data with the shift determined by the distribution 
functions (a),(b) and (c) of Fig. 1 is presented in Fig. 2. The reference Fermi level was fixed 
a.t ep = —0.45 eV, and the hole concentration 6 = 1 — f was calculated by assuming that 
/ = -Dp|/^i=/^2=eF+'5/i/-^pUi=Ai2=eF- The Fermi-Dirac distribution produces a linear decrease 
(Fig. 2-a), while ('^)(^) (Fig- 2-b), and (Fig- 2-c), both lead to a dependence A/i ~ 5^. 

The experimental values can be fitted by a quadratic polynomial (Fig. 2-d). 

The suppression of the chemical potential shift for small 6 compared to its value at 
overdoping is taken as an indication for the opening of the pseudogap [Q. The present 
calculations indicate that the shift produced by an entangled Fermi-Dirac distribution can 
be very well simulated by introducing a fictitious gap parameter. The entangled distribution 
leads also to a quadratic dependence of AyU on 5, as seen in experiment, in contrast to the 
standard Fermi-Dirac distribution. However, although overall the results are in a reasonable 
agreement with experimental data, the simple model considered here cannot explain the 
complete suppression of the chemical potential shift for small 6. 



V. SUMMARY AND CONCLUSIONS 



The effective number of microscopic degrees of freedom in a many-particle system may 
change when classical and quantum aspects coexist, and such variations are reflected in the 
expression of thermodynamic and statistical functions. For a classical particle, there are 
few translational and rotational degrees of freedom. The classical motion can, in principle, 
be described in terms of an infinite number of fictitious harmonic oscillators in a Fourier 
transform, but these oscillators have no statistical relevance. The same holds to a certain 
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degree in the case of a quantum object. The energy eigenstates are normahzed, and have an 
oscillatory time evolution. However, as shown by the quasiclassical picture of the old quan- 
tum mechanics, these eigenstate oscillators provide only a partition of the low-dimensional 
classical phase space into cells of the size of the Planck constant h rather than representing 
new physical degrees of freedom. However, new degrees of freedom are represented by the 
occupation numbers of these cells. 

The quasiclassical dynamics of the occupation numbers in a system of fermions at fi- 
nite temperature was derived here using Grassmann coherent states. We have shown that 
these states admit a real parametrization in which they represent su{2) coherent states 
(Appendix). With respect to this parametrization the expectation value of the occupation 
number obeys Hamiltonian dynamics (Sect. II), and the Grassmann coherent states satisfy a 
closure relationship (Appendix). This relationship allows us to recover the standard Fermi- 
Dirac distribution function (Sect. II) and to calculate the particle density in entangled Fermi 
systems at finite temperature (Sect. HI). 

It is interesting to note that within the thermo-field dynamics formalism ||29|, the trace is 
calculated as a "thermal vacuum" expectation value. This "temperature-dependent vacuum" 
is defined in terms of the composite system consisting of the original physical system and 
an identical but fictitious quantum dynamical system. When the fermion operators of the 
fictitious system are replaced by Grassmann variables, the "thermal vacuum" takes the form 
of a Grassmann coherent state. 

The entangled Fermi-Dirac distribution can be relevant in the description of coupled 
2D-electron systems. As an example, in Sect. IV we have considered the case of two 
adjacent CuO planes in the high-T^ superconductor La2_a;Sra;Cu04. We have shown that 
entanglement can produce significant deviations from the standard Fermi-Dirac distribution 
in the normal state. It also leads to a shift of the chemical potential with the concentration 
that is very close to one of a superconducting distribution at finite temperature. The results 
are in reasonable agreement with experimental data, indicating that a more detailed study 
may be important in understanding the pseudogap-related phenomena. 
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Appendix 



|0) . (39) 



The Grassmann coherent states of Eq. (|^) can also be expressed by the action of a unitary 
operator on the vacuum as 

|r) = exp ^(rfc4-afer*; 

The relationship between the Grassmann variables and (k can be obtained by observing 
that for each k, the operators 

Tk T* 1 

S+ = 1 — f4 , S_ = flfcT"^ , 5'o = alttk - - , (40) 

Tfc \Tk\ 2 



with |rfc| = y/r^ ■ generate an su{2) algebra, and that 

exp (rfc4 - a^r*) |0) = exp [\n\ (3+ - 5_)] |0) (41) 

/ tan \ tA t\ I \ 
= cos |rfc| exp — — ^r^aj, |0) . 



Thus, the expression ( |39D has the form of the coherent state (^ with (k = Tfc|Tfc|^^tan \Tk\. 

With respect to the real parametrization introduced in Sect. II, the Grassmann coherent 
states provide a resolution of the identity operator X of the many-particle Hilbert space such 
that X = cpV with 

v = Jdn 10 (CI , (42) 

where cp is the normalization constant and dQ is a volume element on the manifold of 
Grassmann coherent states. 

If the model space contains only one fermion state {L = 1), the Fock space consists 
of the vacuum |0) and a one-particle state |1) = a| |0) , and \() = (|0) — Ci = 
V^(|0) + |1) Ci) with Ci = ziwi, zi = pie''^\ and = 1/(1 + pj). In this case 

{ni\0 = VliCi\ (43) 

where rii = 0, 1. The variables ^i, 0i are related to pi and qi defined in Eq. (|lT]) by a 
canonical transformation, and if the volume element is chosen as dQ = dpidqi/h = d^idipi, 
then the off-diagonal terms are eliminated by the integration over 0i, leaving 

(z| V \k) = 2t,5,u j dUii^io + Plw^wl5a) (44) 
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Because pi = — 1 and 



r«i= r«i(^-i)=^ (45) 

Jo Jo ^1 ^ 



then {i\ V \k) = nSik, which shows that cp = I/tt. 

Similarly, it can be shown that in general, the overlap between the Grassmann coherent 
state of Eq. (0) and a many-body state \[n]) = \nin2...nL) = (al)*^! (a2)"^...(a^)"^ |0), 
rii, n2, ...riL = 0, 1, is, up to a ±1 factor, 

L 

{nin2...nL\0 = l[VlkC- (46) 

k=l 

By choosing d^fl = Y[k=i dikd4>k, we get 

{nin2...nL\ j d^Vt\Q {C,\mim2...rni) = 'n^5[n][rn] , (47) 

which proves the decomposition of the identity 

J=^| d^n\0{C\ . (48) 

Using this formula we can express the trace of a many-body operator A as an integral over 
Grassmann coherent states. Thus, 

Tr(A) = Tr(AJ) = {[n]\AI\[n]) (49) 
= d^(]S[„](M|A|C) (C|[n]). 

However, because the \ [n\) constitute a complete set of many-particle configurations in the 
selected subspace, Yli[n] (CIW) {V^]\ = (Cl^ and 

Tr(A) = -l I rf^r](C|A|C). (50) 
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Figure captions 



Fig. 1. Planar electron density as a function of energy, (a) normal Fermi- Dirac distribution; 
(b) superconducting at finite temperature; (c) entangled Fermi-Dirac. 

Fig. 2. Chemical potential shift as a function of hole concentration. (*) experimental 
values; (a) normal Fermi-Dirac distribution; (b) superconducting at finite temperature; (c) 
entangled Fermi-Dirac; (d) quadratic polynomial fit of the experimenal data. 
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